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Expanding  che  forced  Rayleigh-Plesset  equation  and  its  initial  conditions 
to  the  second  order  in  e,  we  find  that  the  zeroth-order  problem  is  the 
well-known  autonomous  nonlinear  equation  with  nonhomogeneous  initial 
conditions,  giving  free  oscillations  of  a  typical  nucleus.  This  first- 
order  system  is  a  nonautonomous  linear  system  with  homogeneous  initial 
conditions  which  governs  the  forced  bubble  growth.  The  second- order 
system  consists  of  a  linear  autonomous  differential  equation  and 
homogeneous  initial  conditions.  It  is  needed  to  establish  integrab il ity 
conditions  for  the  first-order  solution. 

The  first-order  solution  is  left  for  future  research  and  the  zeroth-order 
problem  is  analyzed  in  the  phase  plane.  Then  a  novel  approximate 
integration,  r  =■  r(u),  is  given  in  terms  of  elliptic  integrals  and 
functions.  We  were  not  able  to  invert  this  solution  and  so  the  inverse 
u  «  u(r)  is  found  numerically.  These  data  are  then  used  to  find  an 
analytical  approximation  for  use  in  future  first-order  calculations. 
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Nomenclature 

Cp  pressure  coefficient  -  (p  -  po)/(l/2)pVo2. 
Cp3  pressure  coefficient  =  (ps  -  po)/(i/2)pVo2 
D  hemispherical  headform  diameter. 

1  air  content  parameter  »  pa/(2a/R0). 

H  maximum  laminar  separation  bubble 
height. 

K  cavitation  number  -  (po  -  pv)/(l/2)pVQ2 

p  static  pressure  at  any  point  in  the  flow. 
pa  saturation  dissolved  air  pressure. 

p0  free-stream  static  pressure. 

ps  static  pressure  in  laminar  separation 
bubble 

pv  vapor  pressure  of  water. 

Q  from  bubble  radius  normalizes  (1  ♦  Q), 
where  Q  ■  (Wr)K/8. 

R(t)  bubble  radius  at  any  instant. 

Re  diametral  Reynolds  number  =  VqD/i/. 

Rq  free-stream  nucleus  radius. 

r  dimensionless  radius  =  R(t)/R 

O 

p  density  of  water 
s  arc  length  along  headform  meridian. 

0  coefficient  of  surface  tension, 
t  "laboratory  time." 
r  "bubble  time  “ 

Tn  slow-time  scales  =  ent,  n  =  0,  1,  2,  ... 


u  normalized  bubble  radius  =  r/(l  *  Q) 

V  free-stream  cutoff  speed. 

CO  r 

VQ  free-stream  velocity, 
i /  kinematic  viscosity  of  water. 

We  diametral  Weber  number  -  VQ/ -/ o/(pbV 
Wr  radial  Weber  number  =  V  N o/(pRQ). 

Introduction 

This  paper  reports  some  new  beginnings 
in  our  efforts  to  apply  the  Rayleigh-Plesset1 
theory  of  cavitation  bubble  dynamics  to  the 
prediction  of  cavitation  inception  on  sub¬ 
merged  bodies  in  flowing  water  at  ordinary 
temperatures.  For  cavitation  initiation  under 
these  conditions,  we  suppose  that  events 
leading  to  inception  are  sufficiently  infrequent 
so  that  one  can  model  the  dynamics  of 
cavitation-bubble  growth  as  though  each  bub¬ 
ble  were  a  small  isolated  isothermal  spherical 
bubble  containing  a  mixture  of  air  and  water 
vapor 

Thus  in  the  interests  of  computational 
simplicity,  we  neglect  the  presence  of  neigh¬ 
boring  walls  and  of  adjacent  bubbles  We  also 
assume  that  cavitation  is  initiated  from 
invisible  air  microbubbles  of  radius  R0  which 

are  distributed  throughout  the  flew  Past 
measurements2  have  shown  that  the  pop¬ 
ulation  of  such  cavitation  nuclei  is  highest  for 
bubbles  having  a  radius  of  only  a  few 
microns  and  that  their  numbers  rapidly 
become  smaller  as  their  sizes  increase.  Of 
course,  we  know  from  potential  theory  that 
the  minimum  pressure  in  a  ,.ow  occurs  on 
the  boundary  and  it  is  in  just  this  region 
that  inception  is  observed.  Nevertheless  we 
shall  account  for  the  presence  of  the 
submerged  body  upon  the  bubble  dynamics 
chiefly  by  its  pressure  distribution  and  its 
boundary  layer,  retaining  the  simplification 
that  the  bubble  motions  show  no  other 
influence  due  to  the  wall.  Even  then  we  shall 
not  be  reticent  about  introducing  further 

•Superscribed  numbers  indicate  citations  in  the 

References 


/ 


30  July  1986 
BRPiBBB:  Ihz 


advantageous  approximations  for  these  and 
related  flow  or  bubble  properties. 

Moreover,  we  shall  consider  only  cavitat- 
mg  flows  which  fairly  recently  have  received 
additional  intensive  experimental  investiga¬ 
tion3,4'5.  Thus,  we  are  relieved  of  the  consid¬ 
erable  effort  required  to  calculate  the  noncav- 
itating  flow  properties  and  we  will  also  have 
a  body  of  laboratory  data  which  can  be  used 
to  test  the  findings  of  the  bubble  dynamical 
calculations 

Past  attempts  to  use  "simple"  bubble  dy¬ 
namics  for  the  prediction  of  cavitation 
inception  and  the  related  scaling  laws  have 
generally  met  with  disappointment.  Therefore 
it  is  useful  to  confine  our  attention  to  our 
belief  that  the  chief  difficulty  involves  the 
proper  use  of  bubble  dynamics.  By  placing  our 
trust  in  reliable  experimental  data  about 
other  things,  we  need  not  be  distracted  by 
the  need  to  predict  every  single  aspect  of  the 
flow  from  first  principles;  although  we  look 
forward  to  the  day  when  this  achievement 
will  be  realized  with  the  help  of  computa¬ 
tional  fluid  dynamics. 

From  the  literature,  one  finds  that  the 
most  extensive  and  consistent  set  of  measure¬ 
ments  has  been  made  by  a  number  of 
investigate,  on  axially  symmetrical  head- 
forms.  Of  these  recent  measurements,  those 
data  pertaining  to  cavitating  and  noncavi- 
tating  flows  around  hemispherical  headforms 
seem  to  be  the  most  thoroughly  reported. 
Consequently,  we  have  chosen  that  flow 
configuration  as  the  basis  for  the  present 
analysis  :n  spite  of  the  fact  that  in  the 
laboratory  this  flow  generally  exhibits  a  short 
laminar  separation  bubble  downstream  of  the 
minimum  pressure  point  as  illustrated  sche- 
maticaily  ;n  Fig  I.  This  feature  offers  some 
advantages  and  some  computational  com¬ 
plexities  in  the  present  instance,  and  it  suffers 
from  the  fact  that  the  kind  of  cavitation 
inception  which  is  associated  with  laminar 
separation  is  not  found  at  Revnoids  numbers 
usually  encountered  on  ship  propellers  and 
other  fuil-scaie  submerged  bodies  Evidently 
we  are  dealing  with  a  form  of  cavitation 
which  is  a  laboratory  artifact  Even  so.  the 
experimental  data  pertaining  to  the  so-caiied 
bubble-ring  cavitation  are  of  a  quality  and 
consistencv  which  permits  them  to  be  most 


useful  to  test  the  efficacy  of  various  theo¬ 
retical  bubble-dynamical  formulations 


Figure  1.-  Laminar  separation  on  a 
hemispherical  headform. 


Theoretical  Background 

In  a  previous  analysis  of  bubble  ring  cavi¬ 
tation6,  the  senior  author  has  suggested  a  se¬ 
quence  of  key  events  that  could  lead  to  the 
inception  of  bubble-ring  cavitation.  The  en¬ 
couraging  agreement  between  observed  and 
calculated  inception  cavitation  numbers  sug¬ 
gests  that  those  basic  findings  which  result 
from  a  consideration  of  experimental  obser¬ 
vations,  suggests  that  the  following  postulates 
could  be  true,  at  least  in  their  major  respects. 

(a)  Thus,  it  is  supposed  that  a  typical  free- 
stream  nucleus  enters  the  boundary  layer  at 
some  point  ahead  of  the  minimum  pressure 
on  the  body  and  moves  with  a  convective 
speed  characteristic  of  that  flow  region. 

(b)  After  arriving  at  tb-  point  on  the  body 
where  the  absolute  value  of  the  pressure 
coefficient  equals  the  cavitation  number,  the 
nucleus  undergoes  vaporous  growth,  reaching 
a  maximum  radius  down  stream  of  the 
minimum  pressure  point 

(c)  The  /apor  bubble  then  becomes  fixed  at  a 
point  inside  the  lam.nar  separation  bubble 
ahead  of  the  turbulent  reattachment  region 
It  can  then  experience  further  growth  by  air 
diffusion  from  the  water  into  the  cavitation 
bubble  because  the  static  pressure  in  the 
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laminar  separation  bubble  is  generally  greater 
than  the  vapor  pressure  and  further  vapor¬ 
ous  growth  wii!  not  occur 

(d)  This  gaseous  growth  ends  when  the  cavi¬ 
tation  bubble  radius  becomes  large  enough  to 
interact  with  the  free-shear  layer  at  the  edge 
of  the  laminar  bubble. 

(e)  This  interaction  causes  the  cavitation 
bubble  to  be  removed  from  the  separation  re¬ 
gion  and  to  move  with  the  local  flow  into  the 
turbulent  reattachment  region  at  the  end  of 
the  laminar  bubble.  Turbulent  eddies  in  the 
reattachment  region  produce  additional  low 
pressures  in  the  water  which  can  promote 
further  vaporous  growth. 

(f)  In  any  case,  the  intense  shearing  action 
within  this  turbulent  flow  may  distort  and 
tear  the  cavitation  bubbles  and  this  leads  to  a 
small  white  band  or  cloud  of  small  bubbles 
which  is  easily  seen  and  heard. 

This  visible  state  of  inception  or  desinence 
is  called  bubble-ring  cavitation.  It  is  not  to  be 
confused  with  a  more  developed  cavitation 
state,  in  which  the  laminar  separation  region 
is  filled  with  little  attached  glassy  cavities 
This  latter  form  is  actually  a  miniature  free- 
streamline  or  cavity  flow  and  it  is  called 
band  cavitation 

Once  the  gaseous  growth  in  the  laminar 
separation  zone  is  established,  the  conditions 
for  the  onset  of  bubble-ring  cavitation  are 
present  and  all  events  experienced  by  a 
typical  nucleus  in  its  later  stages  of  develop¬ 
ment  occur  automatically.  The  separation 
zone  does  insure  a  favorable  environment  for 
gaseous  growth  because  the  water  is  definite¬ 
ly  supersaturated  in  this  region.  It  is  also 
found  that  the  most  favorable  condition  for 
gaseous  growth  occurs  when  the  vapor  bubble 
just  attains  its  maximum  radius  and  simul¬ 
taneously  arrives  at  the  laminar  separation 
zone.  The  accelerations  in  the  liquid  surround¬ 
ing  the  cavitation  bubble  increase  the  degree 
of  supersaturation  at  the  bubble  wall  over 
that  which  would  exist  if  the  accelerations 
were  negligible.  Therefore  the  probability  that 
the  cavitation  bubble  can  become  stabilized  at 
a  fixed  point  in  the  laminar  separation  zone 
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would  be  greatest  at  this  instant. 

Consequently,  the  process  of  initial  va¬ 
porous  growth  is  the  only  process  of  bubble- 
ring  cavitation  needing  detailed  analysis.  A 
legitimate  question  to  be  asked  at  this  point  is 
why  not  let  a  computer  do  the  analysis  for 
us  by  means  of  a  standard  numerical  method 
such  as  a  Runge-Kutta  routine7  After  ail,  this 
problem  involves  only  an  ordinary  pulsed 
Rayleigh-Plesset  differential  equation. 

There  are  two  important  reasons  for 
embarking  upon  a  program  of  analysis.  First, 
cavitation  inception  is  characterized  by  critical 
values  of  the  flow  parameters  such  as  the 
free-stream  velocity,  VQ,  the  headform  diam¬ 
eter,  D,  the  cavitation  number,  fC  =  (p0  -  pv)/ 
(i/2)pVQ2,  the  free  stream  static  pressure,  p0, 
the  liquid  vapor  pressure,  pv,  density,  p,  vis¬ 
cosity,  p,  and  surface  tension,  o,  as  well  as  the 
dissolved  air  content  and  also  perhaps  by  a 
typical  nucleus  radius,  R0,  or  a  probability 

distribution  of  nuclei  sizes. 

It  is  our  aim  to  find  critical  values  of 
these  parameters,  singly  or  combined  as  in 
the  cavitation  number,  Reynolds  number 
and/or  Weber  number,  etc.,  which  define 
cavitation  onset.  And  of  course,  the  form  of 
the  forcing  function  in  any  particular  instance 
will  also  play  a  role  in  all  of  this.  When  one 
considers  the  large  number  of  these  para¬ 
meters,  it  is  apparent  that  successive  numer¬ 
ical  integrations  for  the  many  cases  of  inter¬ 
est  is  not  likely  to  be  fruitful  from  the  view¬ 
point  of  basic  understanding  of  any  inception 
process  even  if  the  work  were  to  be  done. 

Second,  the  nonlinearity  of  the  Rayleigh- 
Plesset  differential  operator  makes  the  task  of 
relating  critical  flow  parameters  to  the  calcu¬ 
lated  growth  curves,  R(t),  very  demanding 
because  of  the  extreme  accuracy  required6. 
For  example,  it  has  been  known  for  some 
time  that  this  differential  equation  can  have 
solutions  of  different  classes  or  forms, 
depending  upon  the  values  chosen  for  the 
basic  pulse  forcing-function  parameters.  Figure 
2  illustrates  the  kind  of  things  that  can  occur 
This  figure,  in  two  parts,  shows  transient  re¬ 
sponses  of  an  isothermal  nucleus  to  parabolic 
pulses  of  various  pulse  heights,  p,  as  indicated, 
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but  having  a  fixed  duration  9  =  0.4."  It  is  seen 
that  illustrations  (a)  and  (b)  in  Fig.  2  have  dif¬ 
ferent  scales  on  both  axes  and  that  the 
curves  in  (b)  start  at  values  of  time,  r,  where 
the  curves  of  (a)  stop.  The  different  scales  of 
the  two  drawings  prevent  us  from  joining 
them  into  a  single  illustration.  The  reason  for 
this  scale  difference  can  be  appreciated  from 
the  fact  that  the  ordinate  of  each  drawing  is 
the  natural  logarithm  of  the  dimensionless 
radius,  r  =  R(t)/RQ.  As  we  look  at  some  of 

these  curves,  it  is  clear  that  r(t)  varies  by  as 
much  as  a  thousand  fold. 
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Figure  2  -  Isothermal  cavitation  bubble 
growth  curves  in  response  to  single 
paraboiic-puise  forcing  functions 

(a)  Shorter  growth  periods. 

lb)  Longer  growth  periods 


‘This  notation  should  not  be  confused  with  laler 
usa^e  in  which  some  of  these  symbols  will  have 
different  meanings  These  unpublished  results  were 
obtained  in  about  1954  by  the  senior  author  from  a 
series  of  Reeves  analogue  computer  runs  at  JPL. 
They  were  obtained  under  the  direction  of  M  S 
Plessel  The  work,  was  sponsored  by  the  late  Phillip 
Eisenburg,  then  head  of  the  Mechanics  Branch, 
ONR. 


On  the  other  hand,  some  cup/es  show 
essentially  no  net  growth.  The  bubbles  in  this 
case  experience  a  nonlinear  oscillation  having 
a  total  displacement  of  about  four  as  can  be 
seen  in  (b)  for  p  =  0.16  If  the  vaiue  of  p  is 
changed  from  0  16  to  0  20,  the  vapor  bubble 
shows  explosive  growth.  Near  t  =  1000  the 
growth  curve  seems  to  approach  an  asymp¬ 
totic  growth  rate  having  a  nearly  linear  r(t). 
Evidently,  these  contrasting  solutions  give  us 
a  glimpse  of  the  kinds  of  responses  that  can 
occur  and  we  are  left  to  wonder  about  the 
precise  value  of  p,  between  0.16  and  0.20, 
which  marks  the  transition  between  these 
two  radically  different  solutions. 

This  situation  would  not  be  particularly 
bad  if  it  were  not  for  the  fact,  already  noted, 
that  in  some  cases  very  small  changes  in  the 
various  parameter  values  can  have  extraor¬ 
dinarily  large  and  unexpected  effects  on  the 
nature  of  the  bubble  growth6.  Since  the  esti¬ 
mation  of  these  critical  parameters  requires  a 
certain  amount  of  a  priori  knowledge,  we 
conclude  that  even  a  crude  analysis,  if  it  can 
be  properly  executed,  can  be  better  than  none 
at  all.  At  this  time,  all  we  can  offer  regarding 
the  appropriateness  of  one  or  another  analyti¬ 
cal  approach  is  how  well  it  compares  with 
experiment. 

Recently,  one  such  analysis  for  bubble¬ 
ring  cavitation6  has  shown  rather  good 
agreement  between  analytical  estimates  and 
observed  trends  of  cavitation  inception  (  or 
desinence,  in  the  experiments  ).  These  results 
are  shown  in  Fig.  3  where  the  incipient 
cavitation  number  is  plotted  for  a  range  of 
free-stream  velocities.  The  figure  shows  the 
calculated  values  to  be  sensitive  to  the  ratio 
of  maximum  vapor  bubble  diameter,  2R|max  , 

and  to  laminar  separation  bubble  height,  H 
Better  agreement  between  calculated  and  ex¬ 
perimental  findings  could  be  shown  had  we 
permitted  the  ratio,  9  =  2R|  aX/H,  to  exceed 

unity  But  since  the  theory  has  been  greatly 
simplified,  it  seems  that  agreement  between 
experimental  and  theoretical  trends  is  the 
main  point  to  be  exhibited.  The  basic  infor¬ 
mation  on  laminar  bubble  heights  used  here 
is  a  synthesis  of  available  experimental  data 
as  shown  in  Fig  4 

The  basic  scaling  equation  seen  in  Fig.  3 
was  found  to  be 
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Figure  4.-  Laminar  separation  bubble  heights 
on  a  hemispherical  headform. 

K  *  Cps-  (Re°-79/13.8750We2) 
x(  1  -  i(R0/55.59D)2  Re1-58  ]  ,  (  1  ) 

where  the  Reynolds  number.  Re=V0D/V,  is 
based  on  headform  diameter,  D,  and  the 
kinematic  viscosity  is  u=ti/p.  The  Weber  num¬ 
ber,  We  =  VQ//i a/( pD)  is  also  based  on  head- 
form  diameter.  The  dissolved  air  content  pa¬ 
rameter,  7  =  ps/(2o/R0),  is  seen  to  be  the  ratio 

of  the  free-stream  saturation  value  of  the 
partial  pressure  of  dissolved  air,  p  to  the 
free-stream  nucleus  surface-tension  pressure, 
2a/R0.  In  most  laboratories  this  partial  pres¬ 
sure  in  the  free  stream  is  measured  by  a  Van 
Slyke  apparatus.  The  quantity  Cpj  is  the 

pressure  coefficient  on  the  headform  at  the 
point  of  laminar  separation 
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An  interesting  further  aspect  of  bubble¬ 
ring  cavitation  following  from  Eq  (  1  )  is  that 
the  dissolved  air  content  has  practically  no 
influence  on  the  inception  cavitation  number 
when  the  temperature  is  held  constant.  This 
same  conclusion  had  already  been  reached  by 
Carroll  and  Holl4  as  a  result  of  their  experi¬ 
ments.  The  reason  for  this  situation  indicated 
by  Eq.  (  1  )  can  be  seen  from  the  fact  that  in 
the  square  brackets  the  quantity  RQ2  is 

present.  This  quantity  is  very  small  compared 
to  unity  and  so  it  is  found  that  the  entire 
term,  ,f(Ro/55.50D)2Re!  58,  is  much  less  than 

unity  over  the  Reynolds  number  range  of 
interest  here.  Consequently,  one  makes  very 
little  error  by  writing 

K  -  -C  -  (Re°-79/13.8750We2)  (la) 

Ps 

in  place  of  Eq.  (  1  ). 

We  can  not  conclude  however,  that  the 
theory  would  still  be  useful  if  we  were  to  put 
1=0.  If  the  water  contains  no  dissolved  air, 
postulates  (d)  through  (f)  would  have  no 
physical  basis.  Then  there  can  be  no  gaseous 
bubble  growth  and  it  is  doubtful  that  bubble 
-ring  cavitation  would  exist.  Thus,  there  must 
be  some  air  in  the  water  if  bubble-ring 
cavitation  is  to  be  expected. 

We  have  observed  that  a  typical  free 
-stream  nucleus,  which  participates  in  the 
events  leading  to  bubble-ring  cavitation,  has  a 
very  small  size.  That  this  might  be  the  case 
was  inferred  by  matching  the  findings  of  the 
theory  with  a  phenomenon  peculiar  to  bubble 
-ring  cavitation  which  we  will  call  cavitation 
cutoff . 

Cavitation  cutoff  was  observed  first  by 
Carroll  and  Holl  in  their  experiments  and 
found  independently  by  Kodama*  at  about  the 
same  time.  It  was  found  by  all  of  them  that 
there  is  a  lower  limiting  free-stream  "cutoff" 
velocity  for  bubble-ring  cavitation  on  hemi¬ 
spherical  headforms.  At  lesser  speeds  bubble¬ 
ring  cavitation  disappears  altogether,  although 
other  forms  such  as  band  cavitation  may  be 
present.  Our  cutoff-speed  matching  calcula¬ 
tions  indicated  that  the  radius  of  the  typical 
participating  nucleus  in  the  Holl-Carroll  ex¬ 
periments  was  about  5-7  microns.  Aside 


‘Private  communication,  1979 
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from  a  report  to  the  sponsor,  this  theoretical 
development  has  not  been  given  wide 
dissemination.  Therefore,  it  seems  appropriate 
to  present  it  here. 

Capitation  Cutoff 

Basic  Ideas 

As  a  first  step  in  calculations  for  cutoff, 
we  offer  an  interpretation  of  the  experi¬ 
mental  findings  which  can  be  expressed  in 
terms  of  the  present  theory.  One  possibility 
may  be  that  because  the  cutoff  speed  is  a 
relatively  low  speed,  the  cavitation  number 
at  this  speed  should  be  less  than  that  seen  at 
higher  speeds.  The  idea  here  is  that  since  the 
dynamic  pressure  is  less  at  a  lower  speed,  the 
free  stream  static  pressure  must  also  be 
reduced  in  order  to  promote  vaporous  bubble 
growth.  Now,  the  bubble-ring  cutoff  speed 
means  that  if  the  free  stream  velocity  is 
lowered  still  farther  no  matter  what  is  done 
with  the  free  stream  static  pressure, 
bubble-ring  cavitation  is  not  observed,  al¬ 
though  other  forms  may  still  appear.  But  as 
the  free-stream  static  pressure  is  lowered  the 
maximum  vapor  bubble  size  initially  stabi¬ 
lized  in  the  laminar  separation  region  will  in¬ 
crease.  Therefore  we  suggest  that  the  vapor¬ 
ous  cavitation  bubble  is  limited  to  a  greatest 
maximum  size.  According  to  postulate  (d)  of 
the  "theoretical  background"  section  this  size 
must  depend  on  the  laminar  free-shear  layer 
height,  H.  That  is,  there  will  be  little  or  no 
gaseous  growth  phase  at  bubble-ring  cutoff 
speed.  Instead  the  vapor  bubble  is  swept  al¬ 
most  immediately  into  the  tubulent  reattach¬ 
ment  zone.  This  line  of  reasoning  suggests 
that  the  limiting  condition  leading  to  bubble- 
ring  cavitation  cutoff  can  bewritten  as 

aw H  < 2 » 

Figure  4  above  illustrates  that  a  satisfac¬ 
tory  representation  of  the  van  der  Meulen 
and  Arakeri  data  on  laminar  bubble  height  is 
given  by 

H/D  =  111/Re0-79.  (  3  ) 

Since  the  maximum  cavitation  bubble  diam¬ 
eter  is  to  equal  0H  and  if  we  introduce  the 
dimensionless  cavitation  bubble  radius, 


we  can  write 

r  =  55  56D/RnRe0-79.  (  4  ) 

m  0 

The  parameter  0  is  at  our  disposal.  Again,  the 
idea  here  is  that  the  greatest  amount  of 
vaporous  growth  consistent  with  the  key 
postulates  of  bubble-ring  cavitation  theory  is 
limited  by  the  height  of  the  laminar  bubble. 

Rpronlmate  Forcing  Function 

Recall  next  that  the  crucial  simplification 
introduced  into  the  crude  theory  under  re¬ 
view  is  that  the  actual  forcing  function  for 
vaporous  growth  is  replaced  by  one  positive 
step  function  starting  with  a  positive  jump  of 
a  =  -(C„  *  K)/C„  (  5  ) 

where  Cpmmis  the  minimum  pressure  coef¬ 
ficient  on  the  headform  and  Cfl  is  the  dimen¬ 
sionless  coefficient, 

C0=  4O/R0pV02.  (  6  ) 

After  a  short  dimensionless  elapsed  "bubble 
time",  t0  ,  measured  in  units  given  by 


the  first  jump  is  followed  by  a  negative  jump 
of  magnitude 

p  =  (K*Cps)/C0,  (8) 

at  t  -  (ra+  t0).  This  "piecewise  autonomous"  ap¬ 
proximation  is  illustrated  in  Fig  5  .  Adopting 
this  approximation,  we  can  write  the  dimen¬ 
sionless  forcing  function,  F  (t),  as, 

[  0,  0i  t  i  ta  , 

Fa(t)  =  a,  V  <  9  > 

■  P>  WriV 


ffri 


Figure  5.-  Forcing  functions  for  the  "piecewise- 
autonomous"  approximation. 
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Formulation  of  the  Theory 

Next  we  write  the  Rayleigh-Plesset  equa¬ 
tion  as  the  first-order  pair, 

v  *  dr/dt  , 

(2r2)"ld(rV)/dr  =  j/r3  -  1/r  *  Fa(t) ,  (105 

with 

v(0)  -  dr(0)/dt  -  0  and  r(0)  -  1  . 

The  air  content  parameter  i  has  been  defined 
in  the  discussion  of  Eq.  (  !  ).  Multiplying  both 
sides  of  Eq.  (10)  by  2r2  and  integrating 
formally  from  r  =  1  to  r  »  rm,  we  find  that 

rm3vm2=  21.(rm)  *  21  ln  rnT  (  rm2  "  1  <  11  > 

where  r  „ 

‘a  r0  rm 

Ia(i'm)  *  J  r2Fa(r)dr  *j  r2Fa(t)dr  +  J  r2Fa(t)dr 
1  ra  r0 

•  (i/3)[(a  ♦  p)(r03-  0  -§(rm3- 1)]. 

Applying  all  of  this  to  Eq.  (10)  ,  we  get 
rm3vm2  3  <2/3)(Q  *  p)ro3  -  (2/3)a  -  (2/3)0rm3 

+  21  ln  rm  ‘  (rm2  "  •  02) 

Next  one  observes  that  solutions  for 
which  the  bubble  just  achieves  its  maximum 
radius,  rm,  after  Fa(t)  first  equals  requires 
that  phase-plane  trajectories  indicating  un¬ 
limited  vaporous  growth  are  not  admissible 
solutions.  Moreover  one  can  argue  that  the 
demarcation  between  these  unwanted  trajec¬ 
tories  and  those  leading  to  bounded  (periodic) 
solutions  is  the  separatrix  through  the  saddle 
point  on  the  phase-plane  r  axis  at  (rc,0)  .  We 
shall  call  the  saddle  point  location  the  critical 
point,  and  we  shall  denote  the  values  of  p 
and  a  at  this  point  by  acand  . 

It  is  known  that  the  bubble  time,  t,  re¬ 
quired  for  the  radius,  r(t),  to  reach  the  saddle 
point  along  the  separatrix  is  infinite.  But  it 
turns  out  that  F(t)  is  such  that  the  time 
available  in  units  of  "fast"  time  from  c»0  to 
t=tm  is  very  long,  even  though  the  corre¬ 
sponding  period  in  laboratory  or  "slow"  time 
is  very  short  because  of  the  conversion  factor 
between  the  two  shown  in  Eq.  (  7  )  This 
means  that  the  actual  time  available,  as 
defined  by  F(t),  will  place  all  actual  trajec¬ 
tories  very  close  to  the  separatrix.  Moreover 
the  values  of  a  and  p  borne  by  these  trajecto¬ 
ries  will  in  fact,  be  almost  indistiguishable 
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from  acand  pc .  Consequently,  we  make  negli¬ 
gible  error  in  most  situations  of  interest  if  we 
replace  the  actual  a  and  p  by  acand  p  This 

implies  that  we  give  up  the  possibility  of 
matching  the  time  available  for  a  prescribed 
growth,  r  with  the  time  actually  needed  to 

obtain  it.  This  simplification  contains  an  ex¬ 
tremely  small  error  and  it  offers  a  straight¬ 
forward  way  for  us  to  find  r  =  r  ,  a  and  p  . 

7  c  m  c  rc 

In  the  phase  plane  the  saddle  point  is 
found  from  the  conditions  that  vc=Q  and 

[dv/drm](r=rm)  =  0.  From  an  application  of 
these  to  Eq.  (12),  we  find 

0  »  (2/3)(ac  -  pc)r03  -  (2/3)ac  -  (2/3)pcrm3 

+  21  ,n  rrn  -  (rm2  -  1}  (13) 

ana 

pc  -  KVrm2)  -  l]/rm.  (14) 

Equations  (13)  and  (14)  start  us  toward  the 
determination  of  the  cavitation  cutoff  speed. 
They  do,  however,  contain  the  unknowns  rQ 

and  ac,  because  if  for  now  we  regard  VQ=  Vco 
as  known,  and  rewrite  Eq.  (  4  )  as 

rm-  55.50D°-2VO79/(RO  Vco0-79)  ,  (4a) 

then  Eq.  (4a)  determines  rm  provided  RQ  is 
also  regarded  as  being  known.  But  pc  is  deter¬ 
mined  by  Eq.  (14).  Adding  Eqs.  (  5  )  and  (  8  ), 
and  regarding  a  and  p  as  critical  values,  one 
finds  that 

cc  +P=(C.  •  +CD)/Ca.  (15) 

Again  we  consider  the  factors,  RQ  and  V0, 
appearing  in  Cs  as  the  known  nucleus  radius 
and  cutoff  speed.  Then,  ac  is  also  determined 
and  rQ  can  be  found  from  Eq.  (13)  because  all 

other  parameters  in  Eqs.  (13),  (14),  (4a)  and 
(15)  are  known.  These  equations  define  a 
functional  relationship  between  Vcoand  R0  . 

It  may  seem  that  the  foregoing  analysis 
completes  our  study  of  bubble-ring  cavitation 
cutoff  speed.  This  is  not  quite  true  however, 
because  we  have  imposed  no  condition  which 
insures  that  the  value  of  Vcocorresponds  to 

the  lowest  possible  cavitation  number.  The 
following  arguments  are  needed  to  do  this 
and  thereby  to  satisfy  all  physical  conditions 
indicated  at  the  outset. 
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We  start  with  the  observations  that 
situations  can  arise  in  which  the  duration  of 
the  pulse,  t0,  is  so  long  in  terms  of  bubble 
time  that  the  phase-plane  saddle  point,  r  ,  is 
even  less  than  r0.  This  is  almost  as  though 
we  had  let  t0  and  therefore  r0-«,  thereby 
reducing  the  problem  to  the  familiar  autono¬ 
mous  single  step-function  problem  .  In  such  a 
case  the  limiting  values  of  r  will  be  associ¬ 
ated  with  a  new  critical  value  of  a  which 
will  be  less  than  the  critical  value  of  ac  borne 

by  the  separatrix  found  above  and  which  we 
found  to  depend  on  Arguments  only 

slightly  different  from  those  given  previously 
can  be  used  to  show  that 

ac  »  [(Vrc2)  -  l]/rc, 
with  the  analogue  of  (13)  being 

0  =  (2/3)ac(rc3  -1)  ♦  2i  In  rc  -  (rc2  -  1). 
Eliminating  ac  between  this  pair  of  equa¬ 
tions  we  get  a  single  equation  for  the  deter¬ 
mination  of  r  : 

C 

r  ♦  2K 1  -  3  In  r  )  -3  -  2  [(' V r  2)  -  l]/r  =  0  .  (16) 

w  C  C  C 

Because  ac>  I  (3J,  the  present  value  of  rc  will 

be  less  than  the  previous  value.  Therefore  we 
call  it  the  smaller  of  the  two  critical  radii. 
The  closed  trajectories  to  the  left  of  rc  will 
certainly  lead  to  motions  having  short  periods 
and  very  small  amplitudes.  In  order  to  insure 
that  larger  amounts  of  vaporous  growth  are 
produced,  the  phase-plane  trajectories  starting 
from  the  initial  point  will  follow  paths  which 
are  outside  the  new  separatrix  and  they  will 
bear  larger  values  of  a  than  the  ac  just  found. 

And  since  -C  i  is  fixed,  the  K  values  borne 

Plmin 

by  these  trajectories  will  be  less  than  that 
found  on  the  ac  trajectory.  Moreover,  if  we 
put  r0*  rc  we  will  have  reduced  t0  to  its 

shortest  possibile  duration  and  this  requires 
the  lowest  possible  K  vaiue  to  achieve  the 
desired  vaporous  growth  Therefore  we  can 
modify  (16)  to  read 

rQ*  2K 1  -  3  in  r0  )  -  2  ((-f/r02)  -  l]/rQ=  3  ,  (17) 
again,  because  it  is  found  that  the  true  values 
of  a  are  practically  indistinguishable  from  ac. 
For  values  of  r0  in  the  interval  (1,  r  ),  Eq. 
(17)  will  have  one  root  provided  that  Vco  is 


properly  chosen.  The  fact  that  rQ  depends 
parametrically  on  VCQ  suggests  that  we  should 
rewrite  (17)  as 

G(Vco}  *  V  21(1-3  In  rQ) 

-2[(l/r02)  -  l]/r0-3  =  0,  (17a) 

and  solve  for  G  »  0  by  iteration  on  V  because 

J  CO 

Eqs.  (13),  (14),  (4a),  (15),  and  (17)  now  form  a 
closed  system. 

Cutoff  Calculations 

As  indicated  previously,  this  solution  will 
depend  upon  the  value  of  R0  which  is  taken 
as  being  known.  Then  one  can  use  this  value 
and  the  corresponding  value  of  Vco  to  find  C 

from  Eq.  (6).  Once  this  has  been  done  one  can 
find  the  cavitation  number,  Reynolds  number 
and  all  other  parameters  pertaining  to  this 
particular  flow  condition.  A  fairly  extensive 
series  of  calculations  has  been  carried  out  in 
order  to  test  the  preceding  theory  against  the 
cutoff  experiments  of  Carroll  and  Holl  on  a 
2-inch  (0.0508  m)  diameter  headform.  The 
detail  available  to  us  regarding  Kodama's 
experiments  was  insufficient  for  us  to  use  his 
results  to  test  the  theory  further  although  his 
findings  seem  to  be  consistent  with  the  obser¬ 
vations  of  Holl  and  Carroll5.  At  an  air  con¬ 
tent  of  7.8  ppm  and  80'F  (26.67*C),  Holl  and 
Carroll  observed  bubble-ring  cavitation  at  40 
ft/sec  (12  192  m/sec)  but  no  bubble-ring  cavi¬ 
tation  at  the  next  lower  test  speed  of  30 
ft/sec  (9.144  m/sec),  indicating  a  cutoff  speed 
somewhere  between  these  two  values 

Calculated  cutoff  data  for  V„  on  a  2-inch 

(0.0508  m)  diameter  headform  at  a  temper¬ 
ature  of  80®F  (26  67®C)  and  C  -  -  0.630  are 

Ps 

tabulated  below  to  show  the  corresponding 
radius  RQ  of  a  "typical"  participating  nucleus 

at  two  air  contents,  and  values  of  V„  ex- 
pressed  in  feet  per  second. 

Table  1  -  Participating  Nucleus  Size 


Vco 

Air  Content,  ppm. 

7  5 

15  0 

30 

5.65 

4.28 

40 

4.71 

3.61 

t__R 

o__f 

Vco=  9.144  8.  12.192  m/s 
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In  contrast  to  the  findings  of  Eq.  (  1  ),  Table  I 
shows  at  cavitation  cutoff  that  the  nucleus 
size  is  sensitive  to  the  dissolved  air  content. 
For  a  given  cutoff  speed,  the  value  of  RQ 

decreases  with  increasing  air  content.  At  fixed 
air  content,  the  value  of  R0  decreases  with 
increasing  cutoff  speed.  Comparison  of  these 
tabulated  results  with  the  Holl-Carroll  data 
suggest  that  an  appropriate  diameter  of  the 
typical  participating  nucleus  which  will  fit  the 
calculated  results  to  the  experiments  should 
be  between  9  and  11  microns.  This  is  at  the 
small  end  of  experimentally  determined  dis¬ 
tribution  data  in  water  tunnels  and  suggests 
that  there  can  be  many  nuclei  within  the 
size  range  which  may  permit  the  occurrence 
of  bubble-ring  cavitation.  The  calculations 
from  Eq.  (  1  )  leading  to  the  results  shown  in 
Fig.  3  were  carried  out  for  an  RQ  of  5.18 

microns. 

Holl  and  Carroll  also  discovered  a  temper¬ 
ature  dependence  of  bubble-ring  cutoff  cavi¬ 
tation  number  with  Reynolds  number.  Cal¬ 
culations  which  included  the  influence  of  tem¬ 
perature  on  R0,  seemed  to  exhibit  the  same 
qualitative  trends  as  the  observations  but  the 
calculated  results  did  not  display  the  same 
degree  of  temperature  sensitivity,  Fig. 6.  In 
these  calculations  the  physical  parameters 
which  were  taken  to  vary  with  temperature 
were  the  Henry's  law  constant,  the  surface 
tension,  the  vapor  pressure  and  the  kinematic 
viscosity  of  the  water.  Evidently,  other 
factors  should  be  considered  to  complete  the 
theory 


Figure  6  -  Observed  and  calculated  tempera¬ 
ture  effect  on  cavitation  cutoff 
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Because  of  its  several  successful  predic¬ 
tions  so  far,  the  bubble  dynamical  approach 
to  understanding  this  one  type  of  inception, 
argues  for  its  further  refinement  in  order  to 
remedy  the  above  and  other  shortcomings  in 
the  method.  The  most  obvious  topic  for 
further  study  is  the  development  of  a  more 
realistic  treatment  of  the  bubble  dynamics 
which  might  be  capable  of  lessening  the  ob¬ 
served  deficiencies.  The  following  discussion, 
outlining  some  research  in  progress,  summa¬ 
rizes  our  first  attempt  to  address  this  aspect 
of  the  problem. 

The  Refinement 

The  Flaccid  Bubble  and  Initial  Conditions 

A  minor  refinement  to  the  bubble  dy¬ 
namical  theory  underlying  the  preceding  cut¬ 
off  calculations  is  that  we  can  investigate 
possible  benefits  that  a  more  accurate  repre¬ 
sentation  for  the  initial  conditions  of  vaporous 
growth  might  provide.  Generally  this  refine¬ 
ment  is  honored  by  its  neglect.  Then  one  sim¬ 
ply  starts  with  the  undisturbed  free-stream 
nucleus  and  puts  R(0)  *  R0  and  (dR/dtJjl=0-  0, 

because  in  many  cases  of  interest  the 
solutions  for  unlimited  growth,  as  illustrated 
in  Fig  2b,  are  sought  The  precision  of  the 
initial  conditions  is  then  of  little  consequence. 
This  situation  does  not  apply  here  because, 
now  the  vaporous  growth  is  limited  by  the 
separation  bubble  height. 

On  the  other  hand  the  flow  over  a  body 
produces  static  pressures  on  the  surface  which 
differ  from  the  free-stream  value.  In  particu¬ 
lar,  the  static  pressure  at  the  point  on  the 
body  at  the  start  of  vaporous  growth  is  given 
by  Cp-  -  K,  corresponding  to  a  free-stream 

static  pressure  of  p0  For  an  isothermal  flaccid 
bubble,  we  have  pgR(0)3  -  paR03  from  Boyle's 
law.  This  relationship  represents  a  shorthand 
way  of  saying  that  a  free-stream  nucleus 
enters  the  boundary  layer  and  is  conveyed, 
as  an  air  bubble  without  any  vaporous 
growth,  to  the  point  where  Cp^  -  1C,  and  then 
vaporous  growth  starts 

Next  we  can  calculate  the  force  balance 
across  the  bubble  surface  when  it  is  in  the 
free-stream  and  the  outside  static  pressure  is 
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dq.  Then  p0  =  pv*  pa-  2a/RQ  1/  p(0)  is  the 

pressure  outside  the  bubble  at  the  location 
where  vaporous  growth  starts,  the  force 
balance  is  p (0)  =  pv*  p  -  2a/R(c).  The  gas 
pressures  in  these  equations  can  be  found 
directly  and  these  pressures  substituted  into 
Boyle's  law  above.  In  doing  this  we  employ 
the  various  dimensionless  quantities  used  in 
Eqs.  (  1  )  and  (la)  and  we  introduce  the  new 
parameter,  Wr,  the  Weber  number  based  on 
R0  and  write  Wr  =Vq// o/(pRQ)  After  further 

manipulation  of  Boyle's  law  we  get 

[(C  *  K)WR/4]r3  +  r2  -  1  -  (Wr )U4  »  G.  (18) 

P 

Eq.  (18)  is  the  flaccid  bubble  equation  which 
holds  on  the  surface  of  the  body.  When  Cp=  -K 
and  r  =  r(0),  we  can  solve  the  resulting  qua¬ 
dratic  to  get 

rCO)  =  / 1  *  [(Wr)U41  ~  l  *  (Wr)C/8 

=  1  +  Q.  (18a) 

Eq.  (18a)  gives  the  initial  value  of  r  and  we 
write  Q  s  (Wr)K/8 

Our  next  task  is  the  development  of  an 
expression  for  dr/dr  and  its  evaluation  at 
that  point  on  the  body  where  Cp=  -  K  .  The 
known  quantity  on  the  headform  is  C  (  s  }, 
from  which  the  local  flow  velocity  just 
outside  the  boundary  layer  can  be  found. 
Then,  m  keeping  with  the  simplification  that 
the  boundary  layer  is  a  vortex  sheet  we  put 
the  convective  speed,  V(s),  equal  to  one  half  of 
the  outer  fiow  speed.  In  addition,  if  ds  is  an 
increment  of  dimensionless  arc  length  and  dS 
is  the  same  increment  measured  in  suitable 
length  units  we  can  write  ds  =  dS/( D/2).  Thus 

V(s)  -  dS/dt  -  (Dds/2)/dt 

and 

ds/dt  5  (V0/D)v  11  -  C  (  s  )_i  (19) 

Now  we  can  w  rite 

dr/dc  =  (d  r/dCp,'(dCp/ds)(ds/dt)(dt/dt) 

The  first  factor  in  this  product  can  be  found 
by  implicit  differentiation  of  Eq  (19),  the 
second  from  the  equation  for  C  (  s  ),  the  third 
from  Eq  (19),  and  the  fourth  from  Eq  (  7  ) 
After  a  certain  amount  of  manipulation  and 
the  use  of  Cp=  -  K  in  the  appropriate  places, 

we  find  an  expression  for  (dr/dt)i(c .  0)  .  This 

expression  contains  certain  terms  requiring 
r(0)  from  Eq  (18a)  in  order  to  complete  the 


calculation  of  the  initial  bubble  growth  rate. 
When  appropriate  values  of  all  physical  and 
geometric  quantities  are  used  in  the  growth- 
rate  equation,  (dr/di)i(c .  0)  is  found  to  be  ex¬ 
tremely  small  compared  to  other  quantities  in 
the  theory  Therefore  we  make  practically  no 
error  if  we  take  the  second  initial  condition  to 
be  null: 

(dr/dt)|(c,0)-  0.  (20) 

Nevertheless  as  we  shall  see  later,  it  is  of 
conceptual  value  to  understand  that  actually 
this  initial  growth  rate  is  a  very,  very  smail 
positive  number. 


Forcing  Functions 


As  we  have  already  indicated  from  the 
beginning,  vaporous  bubble  growth  starts 
when  C  (  s  )  =  -  K  and  the  dimensionless  bub¬ 
ble  time  or  fast  time,  t,  is  zero.  Let  the  di¬ 
mensionless  arc  length  on  the  body  be  sQ  at 
that  point.  Its  value  can  be  found  from  the 
value  of  K.  But  this  value  varies  and  we  need 
to  allow  for  that  in  defining  the  forcing  func¬ 
tion  because  both  its  amplitude  and  duration 
are  influenced  by  that  vaiue.  On  the  other 
hand  we  can  select  a  likely  range  for  1C  de¬ 
pending  on  the  flow  conditions  or  from  the 
data  of  Holl  and  Carroll.  Therefore  we  can  use 
their  experimental  pressure  distributions  in 
order  to  define  an  average  pressure  distri¬ 
bution  for  the  free-stream  velocity  range  of 
interest  We  can  pick  some  point,  st  say,  at 
which  Cp(  Sj  )  >  -  1C  and  we  can  then  adapt 
Eq.  (19)  in  order  to  define  a  relationship 
between  s  and  a  quantity  t*  that  is  invariant 
with  respect  to  1C.  Clearly,  the  required 
relationship  is 


t*(s)  =  (D/Vo)7(2o/pRQ3) 


dx//[l  -  Cp(  x  )  ].  (21) 


Equation  (21)  is  used  to  carry  out  numer¬ 
ical  integrations  based  directly  on  individual 
points  from  the  Carroll-Holl  data.  Therefore 
one  selects  Sj  to  be  the  first  datum  point  in  or 

near  the  range  of  interest,  C  (s^  -  -  Kmax, 
and  integrates  numerically  point  by  point 
over  the  rest  of  the  points,  s1  >  s,,  in  the  data 


set  The  result  will  be  a  tabie  of  s;  and  Cp(3j) 
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versus  rv,  where  the  index  i  starts  at  1  and 
ranges  over  all  points  indicated  above.  Then 
one  can  plot  the  points  Cp(  tV)  in  (r*l(-Cp) 

coordinates.  Next  he  can  define  displaced  (t.F) 

coordinates  in  which  the  origin  is  positioned 

at  that  value  of  t*  for  which  C  (  t*()  -  -  K. 

a^d  in  which  F  =  -  (C  ♦  K).  If  the  minimum 

value  of  C  is  C  i  „  then  the  greatest  value 
p  p  imin  ° 

of  F  is  Fmax=  (*  cp  Imin*  The  displacement 
of  the  (t.F)  system  with  respect  to  the  (t* ,  - 
Cp)  coordinates  can  be  found  by  means  of 
parabolic  interpolation  in  the  (tV,  -  C  )  plane 

Indeed,  if  the  interpolation  formula  covers  an 
adequate  range  of  cavitation  numbers,  the 
formula  can  be  written  down  once  and  for  all 
and  arranged  to  give  thn  displacements  K  - 
C  (  t*  )  and  t*(K)  -  t*  for  eac:.  value  of  K. 

p  1  1 

used  in  any  particular  computation.  Thus  we 
will  get  a  new  forcing  fun.tion  every  time  K 
changes  its  value. 

In  any  case  a  table  of  points,  t.  versus  F(, 
results  from  these  calculations  as  shown  in 
Fig  7.  Our  aim  has  been  to  preserve  the  Holl- 
Carroll  data  points  themselves  as  far  into  the 
calculation  as  possible  and  we  shall  write  any 
member  of  the  family  of  forcing  functions 
symbolically  as  F(  t.  ;  K  ).  The  faired  curve 
through  the  points  in  Fig.  7  results  from  a 
sequence  of  parabolic  interpolations  between 
experimental  points.  The  next  step  is  to  fit 
the  points  defining  F(  t,  ;  K  )  with  a  suitable 
continuous  function  designated  by  F(  t  ;  K  ). 
This  step  will  be  deferred  to  a  later  point  in 
the  argument,  however.  Having  arrived  at 
the  initial  conditions  and  the  forcing  function 
for  vaporous  growth,  we  turn  next  to  the  dif¬ 
ferential  equation  governing  the  process  and 
to  the  statement  of  the  mathematical  prob¬ 
lem  to  be  addressed. 

The  Initial-Ualue  Problem 

All  ingredients  for  the  statement  of  the 
initial  value-problem  are  given  by  Eqs.  (10), 
(18a),  (20)  and  the  continuous  form  of  the 
forcing  function,  F(  r  ;  K.  )  as  just  described. 

Rewriting  Eq.  (10)  as  a  single  second  order 
equation  and  assembling  the  remaining  re- 


-1  2 1  *  »  - 1 - - - 1 

>  n  3  >35  3.10  3.15  0.20 

LABORATORY  OR  SLOW  T 1  Mf .  1  =  Cr 

Figure  7  -  Hoil-Carroll  forcing-function  data. 


0.2  - 


DIMENSION  l£5S  TIME,  cr 


Figure  8.-  Forcing-function  family  members 
for  various  cavitation  numbers. 

suits  we  may  state  the  initial-value  prob¬ 
lem  as 

r(d2r/dr)  +  (3/2)dr/dt  =  q/r3-  1/r 

♦  [U/2)pV02/(2a/R0)]F(t;  K) 

r(G)  -  7(1  *  (  KWrJLM]  ")  -  1  -  (Wr)K/8  (22) 

=■1*0. 

(dr/dt)kc-or  °- 

Recall  that  the  air  content  parameter  is  de¬ 
fined  in  remarks  immediately  following  Eq.(l). 
Aside  from  the  initial  condition  and  the  fam¬ 
ily  of  forcing  functions,  Fig.  8,  which  are  spe¬ 
cific  to  this  discussion,  Eqs.  (22)  exhibit  the 
customary  dimensionless  form  of  the  initial 
value  problem  for  the  Rayleigh-Plesset 
equation. 

It  is  convenient  to  normalize  the  initial 
condition  on  the  radius  in  Eqs.  (22)  so  that  it 
reads  u(0)  =  1.  If  we  put  r  =  (1  ♦  Q)u,  this 
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substitution  will  certainly  normalize  the  r(0) 
initial  condition  and  the  initial  value  problem 
will  now  read 
u(d2u/dt2)  ♦  (3/2)(du/dt)2 

=  y[(l  ♦  Q)5u3]  -  1/(1  ♦  Q)3u  *  f(t;  D/(l  *  Q)2, 
u(0)  *  1,  (22a) 

(du/dt^.o,*  0, 

where 

f(t  ;  1C)  ■  (d/2)pV02/(20/R0)]F(r;  D. 

These  normalized  equations  are  the  basts  of 
the  subsequent  analysis. 

Multiple  Scales  Representation 

The  possibility  that  a  multiple  scales 
analysis  might  be  useful  for  the  present 
problem  follows  from  the  fact  that  a  fast 
time,  t,  as  defined  by  Eq.  (7),  giving  the 
transformation  between  bubble  time  and 
laboratory  or  slow  time,  t,  can  be  seen. 
Moreover,  the  foregoing  analysis  shows 
clearly  that  the  scale  of  the  slow  time  is 
defined  by  the  convective  process  which 
converts  the  pressure  distribution  on  the  body 
into  the  bubble-dynamical  forcing  functions  as 
indicated  by  Eqs.  (19)  and  (21).  If  we  write  T 
to  represent  the  laboratory  time  scale,  in 
which  t  is  measured,  it  is  clear  from  Eq.  (21) 
that  T  =  D/Vq.  Similarly,  we  can  denote  the 

time  scale  of  bubble  time  by  [r]  and  write  [r] 
=R07pR0/2d  which,  aside  from  a  numerical 

factor  0(1),  measures  the  period  of  free 
oscillations  of  a  cavitation  nucleus.  A  small 
parameter  e  can  be  formed  from  the  ratio, 
[tJ/T.  Thus, 

e  =  [(V0R0)/D]VpR0/2a.  (23) 

A  very  rough  estimate  indicates  that  e  will 
be  about  1/1000.  Evidently  the  conversion  of 
dimensionless  bubble  time,  t,  to  laboratory 
time,  t,  is  ec  =  t  .  Here  the  laboratory  time  is 
taken  to  be  dimensionless  as  represented  by 
the  integral  in  Eq.  (2!)  or  by  the  ratio  t/T, 
for  example.  To  recapitulate,  we  observe  that 
the  two  time  scales  characterizing  cavitation 
inception  in  the  flow  on  submerged  bodies  are 
t  =  er  and  t. 

Tbr  next  step  in  the  analysis  is  to  repre¬ 
sent  the  dependent  variable,  u(r,e),  as  is  cus¬ 
tomary,  by  an  expansion  in  powers  of  e  and 
with  independent  variables 
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Tn  *  enr,  mO,  1,  2,  ... 

Therefore  we  write 

u(t.e)  =  u0(T0,  Tr  T2,  ...)  *  eu,(T0,  T,,  T2,  ...) 

*  e2u2(T0,Ti,  T2,  ...)♦...  (2d) 

Since  we  shall  carp/  out  the  expansion  to 
order  e2,  we  require  only  the  independent 
variables,  T0,  T(  and  T2  in  subsequent  work. 

The  chain  rule  is  used  to  express  the  deriva¬ 
tives,  du/dr  and  d2u/dt2,  appearing  in  Eqs. 
(22a)  in  order  to  insure  that  the  resulting 
multiple  scales  expansion  is  consistent  in  its 
powers  of  e.  For  example,  the  initial  condition, 
u(0)  =  1,  takes  the  form  uQ(0)  -  1  +  eUj(0)  * 
e2u2(0)  =  0  for  all  e.  Therefore,  uo(0)  =  1  is  the 
only  nonzero  initial  condition  among  the  u 
Slightly  more  involved  reasoning  shows  that 
all  initial  conditions  among  the  dun/dt  are 
null. 

Turning  now  to  the  forcing  function  as 
defined  in  connection  with  Eqs.  (22a),  we  see 
that  f(r  ;  K)  actually  depends  on  the  laborato¬ 
ry  time,  t  =  et  =  Tj,  only.  Consequently,  one 
expects  that  it  should  be  represented  as  a 
first-order  quantity,  f(t  ;  K)  -  eg(r  ;  K),  where 
the  function  g  is  designed  to  provide  a  suit¬ 
able  analytic  approximation  to  the  Holl-Carroll 
experimental  pressure  distribution  as  illus¬ 
trated  in  Fig.  8.  A  likely  choice  for  the  form 
of  this  "soft"  forcing  function  is  deferred  until 
later 

It  remains  to  apply  the  method  of 
expansion  discussed  above  and  to  present  the 
system  of  equations,  representing  the  initial 
value  problems  of  zeroth,  first  and  second 
orders  in  e.  The  derivation  of  these  results  is 
tedious  and  we  will  write  them  down  with¬ 
out  further  ado. 

0(e°): 

u0(d2u0/dt2)  *  (3/2)(du0/dO2=  V((l  *  Q)5u03] 

-  1/(1  ♦  Q)5u0, 

uo(0)  -  1.  (25) 

(duo/dr)Kt  =  o)’  0 

CXe1): 

d2ut/dr2  *  [  (3/u0)  duQ/drl  duj/dt 
+  [(4/u0)  d2ug/dr2  *  [9/(2u02)](duQ/dt)2 
*  2/HI  ♦  Q)u0)3  ]  u,  =  [t/(u0(l  *  Q)2)]  g(t  ;  K),  (26) 
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Uj(C)  *  0, 

(du,/dt)|(t  ,0) =  0. 

0(e2): 

d2u2/at2  *  [  (3/u0)  du0/dt]  du2/dr  ♦  (  [(4/u0)x 
d2u0/dr2  ‘  [9/(2uq2)]  (du0/dt)2  *  2/[(]  *  Q)uqP}u2 
--  3{tu1/[(l*Q)un]2}  g(t  ;  K)  -  H(t). 
u2(0)  *  0, 

(du;/dOi,c=0)-  0,  (27) 

where 

H(t)  =  U[2/[(!  *  Q)3u04]  *  [4u,/u0  ]  d2Uj/dr2 

*  bUj/[(l  *  Q)u02]  d2u0/dt2*  [3/2uQ](du1/dt)2 

*  9(u1/u02!(du0/dtXdu1/dt) 

♦  [9uI2/2u03](du0/dr)2  . 
These  three  sets  of  equations  give  the  mul¬ 
tiple  scales  initial  value  problem  to  the  orders 
indicated  and  they  should  permit  the  solution 
to  the  order  e  to  be  found.  We  can  oniy  pre¬ 
sent  the  solution  uQ  here,  however 

This  zero-order  problem  contains  the  only 
nonlinear  differential  equation  among  the 
three  given  above  Aside  from  the 
normalization  and  the  lack  of  a  forcing 
function,  it  strongly  resembles  the  autono¬ 
mous  equation  dealt  with  in  the  cutoff  cal¬ 
culations  above  As  we  shall  see,  the  above 
rather  small  changes  can  make  a  significant 
difference  in  the  nature  of  the  solutions  that 
we  shall  obtain. 

The  Zero-Order  Solution 

Analytical  Background 

For  convenience  we  shall  drop  the  sub¬ 
script  zero  temporarily,  and  write 
u(d2u/'dt:2)  ♦  (3/2'(du/dt)2  *  i/[(l  *  Q)5u3] 

-  1/(1  *  Q)5u, 

u(0)  *  1,  (27) 

and 

(du/dt'l(t-o)  =  0 

It  is  well  known  that  if  we  put  du/dr  »  v  we 
can  replace  the  second  order  differential  equa¬ 
tion  in  Eqs  (27)  with  a  coupled  pair  of  first- 
order  equations  These  read 
v  *  du/dt 

and 

d(uV)/du  =  [21/(1  *  Q)s]/u  -  u/(l  *  Q)3 
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The  second  of  these  equations  has  the  first 
integral 

u3’/2  =  (2l  Jtn  u)/(l  ♦  Q)5  -  u2/(l  ♦  Q)3 

♦  [a  ♦  (i  ♦  or3], 

in  which  the  quantity,  [A  *  (1  ♦  Q)~3],  is  the 
constant  of  integration.  We  can  rewrite  this 
to  read 

u3v2  =  (21  2n  u)/(l  *  Q)5 

-  (u2  -  1)/(1  *  Q)3*  A  ,  (28) 

and  the  condition,  u(0)  =  1,  is  satisfied  Now, 
the  left-hand  side  of  this  first  integral  is 
proportional  to  the  kinetic  energy  of  the 
bubble  motion.  The  right-hand  side  is  propor¬ 
tional  to  the  potential  energy,  -  V,  of  the 
bubble  and  the  constant  A  permits  one  to 
adjust  the  level  of  V.  On  the  other  hand  if 
the  condition,  (du/dtV( .  0)  =  0,  is  also  satisfied 
we  must  set  A  -  0.  Moreover,  had  we  re¬ 
placed  the  forcing  function,  f(t;  K)/(l  *  Q)2,  in 
Eq.  (22a)  with  a  step  function,  F  /(I  +  Q)2,  we 

would  add  the  term,  2[Fc(u3-  1)]/[3(1  *  Q)2],  to 
the  right-hand  side  of  Eq.  (28)  above. 

Phose-Plone  Rnolysis 

Next  we  could  plot  V(u),  including  the 
added  term,  as  shown  in  Fig.  9a.  From  its 
level  lines  and  Eq.  (28),  including  the  Fc  term, 

we  can  find  the  corresponding  phase-plane 
trajectories  as  shown  in  Fig.  9b  This  figure  is 
a  very  familiar  illustration  of  the  global 
nature  of  the  well  known  autonomous 
Tution  and  it  shows  where  the  vortex  and 
saddle  points  are  to  be  found  in  the  phase 
plane.  These  can  be  located  precisely  with  the 
help  of  Liapunov's  method  or  directly  from 
the  potential  energy  curve  The  limiting 
velocity,  ■/[ 2FC/(1  ♦  Q)2]  ,  shown  in  Fig  9  is 

due  to  the  step-function  forcing  function  and 
we  see  that  it  is  naught  in  the  present  zero- 
order  solution  uecause  F  1  0 

C 

While  the  foregoing  discussion  reviews 
knowledge  regarding  general  features  of  the 
zero-order  solution,  it  remains  for  us  to  nar¬ 
row  our  view  and  to  examine  the  problem  at 
hand.  This  we  do  by  considering  the  basic  so¬ 
lution,  Eq.  (28),  with  A  =  F  =  0.  Then  both  ini- 

tial  conditions  are  satisfied  and  this  particular 
form  of  the  differential  equation  (28)  enables 
us  to  study  the  phase-plane  trajectories  and 
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the  position  of  the  singularities  on  the  u  axis 
as  the  parameter  Q  varies  over  the  range  of 
interest.  Figure  (10)  shows  the  results  of  this 
investigation  for  one  value  of  the  dissolved 
air-content  parameter,  j  =  1.4  for  saturated 
water.  Now  the  saddle  point  has  been  lost 
and  ail  trajectories  originate  at  u  =  1  by 
design. 

The  dotted  curve  in  the  upper  potential- 
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Figure  9  -  Phase  plane  and  level  lines  for  a 
step-forced  isothermal  cavitation  bubble. 

(a)  Potential  curve  and  level  lines. 

(b)  Phase-plane  trajectories. 

energy  plot  corresponds  to  the  case  where  the 
trajectory  and  the  vortex  are  contained  at 
the  origin  in  the  phase  plane.  The  motion  is 
then  defined  by  the  solution  u(t)  =  1  and  v(t)  - 
0.  This  null-motion  finding  overlooks  the  fact 
noted  previously  that  the  initial  velocity  is 
actually  an  extremely  small  positive  number. 
But  since  this  model  of  the  oscillations  has  no 
zero-order  energy  source  the  resulting  oscilla¬ 
tions  are  well  represented  by  the  trivial 
solution  just  given. 

The  important  feature  is  that  this 
"critical"  case  represents  the  condition  of  de¬ 
marcation  between  two  physically  distinct 
motions.  Trajectories  on  the  left  represent 
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Figure  10  -  Phase  plane  and  level  lines  for  the 
zero-order  analysis. 

(a)  Potential  curves  and  level  lines  for  various 
values  of  Q  and  one  dissolved  air  content. 

(b)  Phase-plane  trajectories,  originating  at 
u  ■  1,  for  various  Q  values. 

small  scale  flaccid  air  bubble  oscillations 
whereas  those  on  the  right  are  somewhat 
larger  amplitude  oscillations  involving  vapor¬ 
ous  growth,  although  the  gas  in  the  bubble  is 
also  important  in  the  motion.  The  way  in 
which  the  phase-plane  singularities  move 
along  the  real  axis  and  the  size  of  the  trajec¬ 
tories  as  the  value  of  Q  changes  from  0  0  to 
0.3  is  clearly  illustrated  in  Fig.  10.  Mcreover, 
the  potential  energy  minimum  is  found  at 
uv=  (-/!)/(]  *  Q),  which  defines  analytically  the 
vortex  point  locations.  For  any  air  content  we 
can  write  Qcrll  »  (/l)  -  1.  If  Qcril3  0,  then  (-/j) 
*  1  and  this  indicates  the  lowest  value  of  dis¬ 
solved  air  content  for  positive  values  of  Q.  Re¬ 
ferring  to  Eq.  (18a)  which  defines  Q,  we  see 
that  this  limit  would  correspond  to  VQ=  0  or 

to  K  -  0,  because  it  is  physically  unrealistic 
for  the  surface  tension  to  be  infinite 

Of  course,  all  of  these  free  oscillations  are 
of  high  frequency.  It  is  of  interest  to  find  out 
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how  many  cycles  of  free  oscillation  might  be 
found  during  the  time  of  rise  for  one  of  the 
forcing  functions  of  Fig.  8.  This  rise  time  is 
measured  in  units  of  laboratory  time.  Some 
calculations  have  been  made  to  answer  this 
question  and  these  will  be  discussed  later.  For 
now  we  may  say  that  during  a  typical 
bubble  time  interval  between  F(0;  K)  »  0  and 
F(t;  K)  =  F|  the  number  of  vaporous  oscil¬ 
lations  is  measured  in  multiples  of  ten.  It 
may  be  pure  coincidence,  but  the  amplitudes 
of  oscillation  indicated  here  and  those  illus¬ 
trated  in  Fig.  2  appear  to  be  about  the  same 
and  so  the  frequencies  also  should  be  com¬ 
parable  It  should  be  said  however,  that  the 
results  found  here  show  a  weak  dependence 
on  1C  and  VQ  and  these  results,  being  of  order 

zero  in  e,  are  independent  of  the  forcing 
function.  Therefore,  the  apparent  consistency 
between  the  present  and  the  older  results  is 
very  encouraging.  Those  older  results  also 
suggest  how  the  influence  of  the  forcing 
function  grows  to  produce  a  large  vaporous 
growth  which  overwhelms  the  small  scale 
free  oscillations  considered  here.  Further 
exploration  of  these  additional  aspects  of  the 
process  depend  on  the  first-order  solution 
which  remains  a  topic  for  future  research. 

Zero-Order  Analytical  Solution 

Using  the  fact  that  v  *  du/dr  and  solving 
for  v  in  Eq  (26),  we  can  write  Eq.  (28)  as 
(du/dr)2  «  [(2i  in  u)/(l  *  Q)5 
-  (u2  -  0/(1  *  0)3]/u3 

The  bubble  time,  t,  is  then  given  as  the  quad¬ 
rature, 

U 

t  -|x2dx/'/(2’!x  2nxJ/(l*  Q)5-  (x3-  x)/(l  ♦  ~.  (29) 

Generally,  one  evaluates  this  integral  numeri¬ 
cally.  It  is  possible  to  obtain  a  reasonably  ac¬ 
curate  approximate  solution  by  fitting  the 
function  u  in  u  with  a  cubic  polynomial,  au3* 
bu2*  cu  ♦  d.  The  coefficients  are  a  =  -0.050304, 
b  -  0.579693,  c  -  -  0  001871,  &  d  =  -  0  526519. 
Then  the  radicand  above  also  becomes  a  cubic. 
Ax3*  Bx2*  Cx  *  D,  having  coefficients: 

A  ■  [2ia/d  *  Q)5]  -  [1/(1  *  Q)5],  B  -  2ib/(l  *  Q)5, 

C  *  2k/U  *  Q)5*  (1/(1  *  Q)3].  &  D  -  [2id/(l  ♦  Q)3]. 
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As  a  rule  this  cubic  can  be  factored  with 
the  help  of  Tartaglia's  method.  But  first  the 
radicand  should  be  rewritten  as  A  [  x3*  px2» 
qx  ♦  r  1 .  Then  if  we  put  uv=  V^f/(1*Q)  we  have 
p  =  2buv2/(2auv2-  1),  q  =  (2cuv2*  l)/(2auv2-  1) 
and  r  =  2duv2/(2auv2-  1)  . 

Let  the  three  real  roots  of  this  modified 
cubic  be  xt,  x2,  and  xJ(  and  suppose  that  they 

are  ordered  with  respect  to  the  upper  limit  of 
integration  to  read,  xt  <  x2  s  u  <  x,.  The  usual 

laborious  factorization  can  be  shortened  in 
this  case  by  observing  that  the  rewritten 
cubic  always  has  one  root  equal  to  unity. 
Thus  when  one  divides  the  cubic,  x3*  px2*  qx  * 
r,  by  the  factor,  x  -  1,  the  remainder  turns 
out  to  be  2uv2(a  *  b  *  c  ♦  d)/(2auv2-  1).  Then 

when  the  numerical  values  of  a,b,c  and  d  are 
summed  as  indicated,  one  finds  this  sum 
equal  to  -0.000001  and  he  can  consider  this 
result  to  be  zero  for  all  values  of  u^,  within 

the  accuracy  of  these  four  coefficients.  Con¬ 
sequently,  we  need  only  factor  the  quadratic, 
x2*  x(2auv2-  1  *  2buv2)/(2auv2-  1)  * 

2uv2(a  *  b  *  c)/(2auv2-  1). 

These  factors  are 

[(1  05b8uv2-  1)/2(1  *  0.100608uv2)]  {  1  i  41  }, 
where 

Z  •  1  ♦  [4.212144uv2(I  *  0.100608uv2)/ 
(1.0568uv2-  l)2]  . 

These  rools  were  also  found  numerically  for 
the  specific  example  of  uv=  1.1832,  correspond¬ 
ing  to  j  =  1  4  and  Q  *  0  0  .  The  three  roots  are 
x,  =  -  0.9459,  x2  =  1  0  and  x3  =  1.3662  . 

Comparison  of  these  with  the  range  of  inte¬ 
gration  in  Eq  (29)  justifies  the  above  ordering 
of  the  roots  with  respect  to  the  value  u,  in 
this  case  at  least.  Moreover,  except  for  Xj,  the 

roots  depend  on  uv,  that  is  xt  =  Xj(uv).  Also  uv 

will  be  recognized  as  the  parameter  governing 
the  location  of  the  singularities  in  the  phase 
plane  The  root  at  u  *  x2  *  1.0  is  seen  to  be 
the  initial  point  in  the  phase  plane.  The  de¬ 
pendence  of  the  other  roots  on  uv,  permits 
one  to  distinguish  between  the  families  of  tra¬ 
jectories  to  the  right  or  to  the  left  of  the 
initial  point  in  Fig  10  above.  The  junior 
author  has  investigated  this  last  matter  fully 
in  his  thesis. 
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Consequently  Eq.  (29)  can  be  rewritten 
using  roots  xlt  x2,  and  x3,  as 

t  *  /(l  *  Q)J/(2auv2*  1)  J  x2dx/G(x)  , 
where  1  (29a) 

G(x)  *  VTx  -  Xj)  (x  -  x2)  (Xj-  x ) 

and  the  value  of  a  is  given  above  .  With  the 
help  of  Byrd  and  Friedman7,  Eq.  (29a)  can  be 
expressed  in  terms  of  incomplete  elliptic  inte¬ 
grals  of  the  first,  second  and  third  kinds  and 
a  product  of  Jacobian  elliptic  functions.  We 
are  not  aware  of  any  method  for  inverting 
that  complicated  result  so  that  it  will  have 
the  form  u  *  f(t).  Therefore  we  did  repeated 
numerical  integrations  and  varied  the  quan¬ 
tity  uv=  /y(  1  ♦  Q)  over  the  range  indicated 
in  the  phase  plane  analysis  above  and  selected 
an  appropriate  intermediate  value  of  (1  *  Q)3/2 
with  the  result  illustrated  in  Fig.  11. 

The  relation  shown  in  Fig.  11  was  fit  with 
an  approximate  analytical  result,  which  is 
accurate  at  the  origin  and  at  all  maximum 
radii,  rm  It  is  a  simple  expression  and  in  the 

worst  case  it  offers  about  5%  or  13%  accuracy 
between  these  limits  depending  on  whether 
one  uses  u  or  (u  -  1)  in  the  calculation  of 
relative  error.  This  simple  aproximation  has 
the  form 

u  =  A  ♦  B  cos(y)  ♦  C  cos(24>).  (30) 

The  angle  ip  =  <p(uv,i:)  and  the  coefficients,  A,B 
and  C,  which  also  depend  on  uv  contain  all  of 
the  complexity.  The  degree  of  success  achieved 
with  this  formula  is  shown  in  the  comparison 


Figure  11.-  Numerical  results  for  zero-order 
isothermal  cavitation  bubble  growth  on  a 
hemispherical  headform. 
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in  Fig.  12.  Note  that  the  ordinate  in  Fig.  12 
has  a  suppressed  zero  so  that  the  illustration 
emphasizes  differences  between  this  approxi¬ 
mation,  shown  as  dashed  curves,  and  solid- 
line  numerical  results  from  Eq.  (29a).  Equa¬ 
tion  (30)  is  an  important  first  step  toward 
achieving  the  desired  precise  result  for  the 
zero-order  solution.  It  suggests  that  over  the 
parameter  range  of  interest  here,  a  more  re¬ 
fined  representation  of  u0(t,  uv)  can  now  be 

found.  This  final  result  and  other  details  not 
presented  above  are  given  in  the  junior 
author's  thesis. 
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Figure  12  -  Comparison  of  numerical  results 
from  Fig.  11  with  the  fitting 
function  of  Eq  (30) 

Conclusions: 

Perhaps  the  chief  point  needing  emphasis 
with  regard  to  our  review  of  the  so-called 
asymptotic  bubble-ring  inception  theory  is 
that  its  postulates  do  not  allow  for  nuclei 
recirculation  within  the  laminar  bubble  after 
they  have  been  created  by  the  visible  bubble¬ 
ring  cloudlets  in  the  transition  zone.  That  is 
to  say,  we  are  not  concerned  with  desinence. 
Our  interest  is  in  inception  only  In  this  latter 
situation  there  will  be  no  ring  of  cloudlets  and 
we  must  assume  that  the  only  nuclei  source 
is  from  the  free  stream.  The  apparent  good 
agreement  shown  in  Fig.  3  leaves  us  with  the 
hope  that  the  basic  assumptions  of  the  theory 
are  close  to  the  truth. 

The  ability  of  the  cavitation-cutoff 
theory,  as  outlined  above,  at  least  to  show 
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correct  trends  reinforces  the  hope  that  our 
basic  postulates  are  correct.  The  question  is: 
Will  a  more  refined  calculation  as  initiated 
here  improve  the  ability  of  the  theory  to 
predict  cavitation  cutoff  speeds?  At  present 
we  can  not  say  But  more  comprehensive 
experimental  cutoff  data  would  be  of  great 
value  for  testing  this  aspect  of  the  theory  . 

Turning  to  the  multiple  scale  calculations, 
we  can  say  that  an  approximate  zeroth-order 
solution  is  now  in  hand  It  remains  to  find 
the  first-order  solution  and  to  put  these  two 
parts  of  the  solution  into  a  form  which  is 
ready  for  more  refined  scaling  and  cutoff 
calculations.  We  hope  that  research  aimed  at 
completing  this  work  can  be  carried  out  soon. 
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